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The generation of finite energy packets of X-waves is analysed in normally dispersive cubic media 
by using an X-wave expansion. The 3D nonlinear Schroedinger model is reduced to a ID equation 
with anomalous dispersion. Pulse splitting and beam replenishment as observed in experiments 
with water and Kerr media are explained in terms of a higher order breathing soliton. The results 
presented also hold in periodic media and Bose-condensed gases. 

INTRODUCTION 

X-waves, originally introduced by Lu and Greenleaf in acoustics, 0, can be simply described as a non- 
monochromatic superposition of plane- waves with the same wavevector component along a given direction of propaga- 
tion. The resulting pulsed beam can travel undistorted without neither diffraction or dispersion, even in the absence 
of a nonlinear self-action. Such X packets have been the subject of intense research in several fields, as recently 
reviewed in 0. [J Q 

The first experimental evidence of X-wave formation during nonlinear optical processes was reported in ■ Sub- 
sequently, X-related results have appeared in 0, II H EI], including a theoretical analysis and discussion on the 
generation and existence of this class of spatially and temporally localized waves in nonlinear media E3 • 

The above mentioned experimental and theoretical results have stimulated further investigations on the role of X- 
waves in nonlinear pro cesses, as well as their description in terms of linear waves with a convenient envelope approach, 
HI El El El El Hi El El El El El E3 and their extension toquantized fields. Within the active field of 
optical self-invariant pulsed beams (see among others [25L l2fil I27L l28l Efl Efl Ell|). it was recognized that X-waves 
may play a fundamental role in all the nonlinear processes encompassing spatio-temporal mechanics, well beyond the 
realm of quadratic frequency mixing where they had been originally observed. X-waves have been also predicted in 
periodic media and Bose condensed gases. 

X- or progressive undistorted waves in focusing or defocusing Kerr media can be associated to the instability of 
wide beams propagating in the presence of a refractive index linearly dependent on intensity. 33] This effect has also 
been investigated in quadratically nonlinear media [20I Ei) . 

Very recently light filamentation with femtoseconds pulses in water, [34| was explained by the "nonlinear X-wave 
paradigm". [35| The authors pointed out that these spatio-temporal packets can be highly dynamic, i.e. they are 
continuously generated, undergoing splitting and replenishment with an "average" invariant propagation. Such a 
picture points out the robustness of X-waves with respect to the numerous non-trivial effects that may take place in 
3+1D nonlinear processes. On the contrary, self-trapped light bullets, which do not exist in the absence of a nonlinear 
response, are extremely sensitive to the specific model. X-wave dynamics seems st rong ly connected to pulse splitting 
during self-focusing in normall y dispe rsive media, as originally predicted in 1992 [36L l37|. and later investigated by 
several authors. |38l l39i l4Ct l4ll l42l I43L l44l |45| Experimental investigation of splitting/replenishment has been reported 
in 0. 

The generation of X-waves may also be taken as the basis for the interpretation of the conical emission observed 
during the propagation of powerful ultrashort laser pulses in air, or water. [3 El El Ell 

In this manuscript a preliminary settlement of X-wave nonlinear dynamics is attempted. While the forthcoming 
theory is developed with reference to nonlinear optics, as outlined in |32( it applies to Bose Einstein condensation (i.e. 
to "matter X-waves"), as well. The central idea is that, if these wave- packets can be considered somehow as "modes 
of free space" , with a given direction of propagation z, then their nonlinear evolution is essentially two dimensional, 
involving z and time. Thus, under certain conditions, the problem can be reformulated and strongly simplified by 
adopting the well known approach of guided wave nonlinear optics. 

Once an appropriate X-wave expansion is determined, it is not straightforward to write down "coupled X-wave 
equations" , in analogy with coupled mode equations, for the simple reason that these spatio-temporal beams are not 
normalizable. The first step is then to build finite energy solutions, and then use the resulting superposition of X-waves 
to investigate nonlinear dynamics. Section I and II cover the early stages of this approach. In section III nonlinear 
regimes are addressed via a perturbative expansion. This is legitimated by the fact that progressive undistorted waves 



2 



do exist in the absence of nonlinearity. The analysis clarifies in which sense X-waves are robust, and may constitute 
a "paradigm" for nonlinear 3D+1 dynamics. In section IV I consider the highly nonlinear regime. This is possible 
while limiting to a specific X-wave. If the wave packet is spontaneously generated during a nonlinear process (e.g. 
from a bell shaped pulsed beam) the resulting picture, for example the emergence of breathing solutions, provides 
a relevant insight for interpreting reported numerical and experimental findings. The last sections deal with some 
corollaries: compression and chirp of X-waves, and the relation with the integrable ID nonlinear Schroedinger equation. 
The appearance of an integrable model, enforcing the use of categories such as "breathers" , and its effectiveness in 
interpreting experimental and numerical investigations, points to an intriguing connection between nonlinear X-waves 
and solitons. 



X-WAVES FOR THE 3D+1 SCHROEDINGER EQUATION 

The wave equation (at angular carrier frequency uiq) describing paraxial propagation in normally dispersive media, 
at the lowest order of approximation, can be cast as 

id z A + ik'd T A+^Vl y A-^-dTTA = 0, (1) 

where k — won{ujQ) / c and its derivatives provide the dispersion. In a reference system traveling at the group velocity 
of the medium, t = T — k'z, it is 

id z A+±-Vl v A-^-d tt A = 0. (2) 



For the sake of simplicity hereafter I will consider radially symmetric beams, with r = -*/X 2 + V 2 , . 

Progressive undistorted waves, propagating with inverse differential velocity (3 |50f can be retrieved by looking for 
solutions of the form A = ip(t — f3z, r) exp(ifc z z), thus (r = t — (3z) 

- k z ^ - i(3d T ^ - y #v + ^ffV> = o. (3) 

If tp is written as a superposition of monochromatic Bessel beams, Jo(Vkk"ar) exp(— ilut), with a in frequency units, 
the corresponding spatio-temporal dispersion relation is 

-k z -^+^ 2 = ^f. (4) 

In order to have a continuous spectrum along uj, the left-hand side must be positive, thus ensuring the absence of 
evanescent waves. This can be achieved, in the simplest way, by letting k z = —/3 2 /2k", which gives 

(--|r) 2 = « 2 - (5) 

Eq. © implies the existence of two types of X-waves: 

^ / (r,r,/3)=/ e- l( ^ +Q)T J (VFfcar)/ / (a)da (6) 
Jo 

and 

/>oo 

^(r,r,/3)= / e-^- Q )V (\/FW)j\(a)da, (7) 
Jo 

with the corresponding "spectra" f/(a) and f\(a). They will be denoted as "slash" and "backslash" X-waves, because 
of the shape of their spatio-temporal frequency content, as discussed below (see Fig^l. 

A general linear X-wave solution. |5l| traveling with inverse differential velocity f3, is given by 

A x = e -^ z [Vv (t-0z, r, /3) +iP\{t-/3z, r, (3)] , (8) 
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FIG. 1: Sketch of spatio-temporal spectra of "slash" and "backslash" X- waves. 



which can be rewritten as 



A x =e-^ t+ ^ z [i P/ (t-Pz 1 r) + ip\(t-Pz 1 r)], (9) 



with 

f°° 

<p X (r,r) = e^ laT J (Vl^kar)f x (a)da. (10) 
Jo 

The "X" stands for either / or \, and ip x corresponds to X- waves of the Helmholtz equation. Q 
The spatio-temporal spectrum of Ax(r, t, z) is given by the Fourier-Bessel transform pair 



1 f°° f c 
B[A](k ± ,uj,z) = — \ 

27r J-oo JO 



A(r, t, z)Jo(k±r) exp(iujt)rdrdt 



A(r,t,z)= / / B[A](k±,oj, z)Jo(k±r) exp(—iujt)k±dk±dLLi, 

J-oo JO 

it is centered at the shifted central frequency f3/k", determined by the velocity: 



(11) 



B[A x ] = ^f/{^=)8{w 

kj_ Vfc"fc Vk"k k" 

— fit , )S((j-\ ■ )e v fc ' 



k± P \ i(gr + 4fc)» 

N g 2k Vk"k 



(12) 



and it looks like an X in the angle-frequency plane. The two terms in l|12|) correspond to slash and backslash X-waves, 
as shown in figure Since the pulsed beam travels rigidly (in modulus), an X-shaped spectrum determines the 
X-shape in the (r,t) space (roughly, the far field, i.e. the Fourier-Bessel transform, resembles the near field). 
The following relation holds useful for any i/;x(t — (3z, r, /?): 



= {id z + ±.V% - yWx = ~^x. (13) 
Therefore an X-wave can also be defined as a solution of the equation @ of the type A x — C(z, f3)ip x (t — j3z, r) with 

a formulation that will become handy below. 

Such a solution contains, in general, an infinite energy. This is due to the idealized situation (never available in 
Experiments) of a precisely defined velocity (or inverse differential velocity (3). Any finiteness introduced by the 
experimental setup, such as the spatial extension of the sample, will in general fade the spectrum lineshape around 
the X, determining uncertainty in j3. I will show below that this can be described by a packet of X-waves, with 
velocities around a given value and finite energy. Considering an X-wave with a specific velocity is thus as idealized 
as considering an elementary particle of given momentum. [52j | 

Before proceeding further I emphasize that entire more general classes of X-waves can be generated as discussed in 
[r?| or, for instance, by taking k z = — 1 /2k" — k 2 in ©. However, the case considered (k = 0) suffices to represent 
a wide class of beams, as it will be shown in the next section. 
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X-WAVE EXPANSION AND FINITE ENERGY SOLUTIONS 



The general solution of Eq. (01 can be expressed by the Fourier-Bessel spectrum of the field at z = 0, denoted by 
S(w,k±) = B[A](w,k x ,0): 



OO POO 



A(r,z,t) = / S{Lu,k ± )J (k ± r)e l( - k * z - ut) k ± dk ± dLj, (15) 



with k z = -k\/2k + k"uu 2 /2. 

By a simple variable change, the field can be written as a superposition of slash- X-waves, traveling with different 
velocities. Letting 



(16) 



k ± = ^/kF'c 



gives 



while being 



and 



OO 



A(r,z,t)= e-^ z ^ / (t-(3z,r 1 P)dp, (17) 



/ 

OO 







9/(t-0z,r,0)= X / (a,/3)J (\/Ffc"ar)e- l(Q+ ^ )( *- /3z) da, (18) 







X/{a,(3) ee X/[A(r,t,z = 0)](a,/3) ee k a S{a + f^, Vkk"a). (19) 

An equivalent representation is obtained by backslash- X-waves, replacing the equation for uj in O with uj = —a + 
/3/k", and being 

X\(a,j3)=X\[A(r,t,z = 0)](a,/3) = kaS(-a+ ^,Vkk"a). (20) 



The variable change (|L(t|) corresponds to span the (uj, k±) space by oblique (slash) parallel lines in Ijl5(l . 

Eq. Ijl8|l is a formulation of the "X-wave transform", first introduced in |£3j and indicated ^[^(a,^). Hence the 
spatio-temporal evolution, including diffraction and dispersion, can be represented by one dimensional propagation 
of packets with different velocities. 

The energy of the pulsed beam is 



OO /"OO 



£ = 2n / \A(r,t,z)\ z rdrdt = £p{(3)d(3 (21) 

J J — oo J ~oc 

with 

f°° 47r 2 |X/(a,/3)| 2 , 
Ep{0) = J ' £ — da, (22) 

showing that Bp can be taken as the energy distribution function with respect to the inverse differential velocity (3. 

Summarizing, an arbitrary beam can be expanded in a superposition of X-waves traveling with different velocities. 
Conversely, such a superposition can be used to construct new-classes of physically realizable finite-energy X-waves. 
To this extent, orthogonal X waves -first introduced in |54j for the wave equation- are a fruitful approach (see also 
[13] )■ With reference to two (either slash or backslash) X-wave solutions of (0, denoted by Ax and Bx, with inverse 
differential velocities (3 and f3' and spectra / and g, respectively, the inner- product can be defined as the integral of 
B* x Ax with respect to x, y, t, extended on the whole space: 



<B X \A X >= J J J B *x Ax dxdydt = 6(p - p) J™ ^ 9 ^ ^ da. 



(23) 
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If / — fp an( l 9 = fg by defining [p, q = 0, 1, 2, .. and ip 1 ^ is the generalized Laguerre polynomial) 



= y ^j^4 1) ( 2Aa ) Aae ~ Aa ' ( 24 ) 

with A a parameter with the dimension of time.[l5| it is 

- ^ (25) 

and the orthogonality condition (A g = Bx and A p = Ax) holds: 

<A q (r,t,z,0)\Mr,t,z,l3') >= S pq S(0 - 0'). (26) 
In the following the fundamental slash-X-wave, with spectrum 

/o(o) = — Aaexp(-Aa) ; (27) 

7T 



and spatial profile 



/•OO /IT A 

^ 0) =/ / (a)J (x^Fa)e-- s d a = -^- fcfc „ r2 — , (28) 

7 -A) T [1 - T |^ F ]3/2( S _ zA )2 

will be taken as a prototype of the simplest X-wave with finite power (i.e. converging transverse integral of l^^l 2 , 
see also [13 )■ Its intensity profile can be found, e.g., in [sij . 
When X{a,0) = C(/3)f p (a), it is £ p = \C{0)\ 2 , and 

/OO 
|C(/3)| 2 d/3. (29) 
-OO 

The resulting beam is a finite energy X-wave that spreads according to a prescribed velocity distribution function 
C(f3); this corresponds to the existence of solutions with an arbitrary "depth of focus". 
This kind of wave-packet can be written, with reference to slash X-waves, as 



A= C(0,z)tp / ( - q) (r,t- 0z,0)d0, (30) 

J — OO 

(with obvious notation) while being, as above, 

<£ -&-<>■ 

By introducing the Fourier transform pair of C (/?) 

(32) 

C{0,z) = ^S™ oo c{s)e- i & s ds, 
with s in time units (s can be roughly kept in mind as the on-axis temporal variable t), it is from Ij31|l 



.dc k"d 2 c 

l d- z + Yd^ = - (33) 

Hence 3D linear propagation of an X-wave packet in a normally dispersive medium can be reduced to that of a ID 
pulse with anomalous dispersion. The energy is 



/OO 
\c(s)\ 2 ds. (34) 
-OO 



As discussed in |32j . by expanding X(a, 0) into generalized Laguerre polynomials it is possible to express the general 
solution of (J2J by orthogonal X-waves (either / or \). This has been used in [24j] as an approach to field quantization, 
and applied to quantum optical parametric amplification. 



6 



GENERAL PROPERTIES OF X-WAVE PROPAGATING IN THE NONLINEAR REGIME 



Since X-waves exist even in the absence of a nonlinear susceptibility, it is legitimate to resort to a perturbative 
approach. A basic model for nonlinear optical interactions (under standard approximations) can be cast in the form 

id z A + ^l y A - \d tt A = X Vnl{z, t, r), (35) 

where VpfL(z,t,r) is a nonlinear source, weighed by x- After a straightforward expansion of A in powers of %, the 
relevant equation is 135f) at any order, where the right-hand side stems from the solutions at lower orders. I will show 
below that if xPnl(z, t,r) =V(t — f3z,r), the evolution according to l|35|) always provides a spatio-temporal spectrum 
corresponding to a progressive undistorted wave. Thus, if an X-wave is taken as the solution of the linear model 
(x = 0), the nonlinearity has the role of "dressing" that solution, which still continues to exist. Since this result is 
valid at any order, it furnishes a general picture of the propagation of self-invariant beams in the nonlinear regime. 

For an X-wave propagating in a nonlinear medium, V can be interpreted as a function of its field and of its complex 
conjugate. In the case of harmonic generation, V is some power of the pump beam, traveling at the group velocity of 
the fundamental frequency. [5^| What follows can be viewed as a generalization of what discussed in 0] , with the 
inclusion of second order dispersion and for a wide class of nonlinear processes (such as third and higher harmonics 
generation) . 

In a Kerr medium, with a refractive index n = hq + ri2 I, with / = \A\ 2 the optical intensity, Eq. H35[) becomes 

id,A + ^y%A -'f-duA+'^IA^O. (36) 

For the solution at the lowest order (ji2 = 0), it is possible to take either an X-wave packet Ax around (3 or a wide 
pulsed beam with negligible diffraction (with (3 — 0). Higher orders are obtained in the form (|35|l : at the first the rhs 
is proportional to A 2 X A* x . As a result of the following analysis, the correction to Ax is still a progressive undistorted 
wave traveling at the same velocity. Since this can be applied at any order, lincarly-sclf-invariant beams are very 
robust with respect to the nonlinearity. 

By writing A, as a superposition of slash- X-waves, 



A= / C(a u p 1 ,z)J o (Vkk^a 1 r)e-^ ai+0 ^ k "^ t ~ filZ Ua 1 d/3 u (37) 

JO J -oo 

and inserting into (|35|l . one obtains 



oo Jo 



(*-^ ~ ■^ 71 C)J {Vk^a 1 r)e' l ^+^l k ^"^Wi dp l = V{t - pz, r). (38) 
az 2k" 



Taking the Fourier-Bessel transform of Ea. (|38|l . multiplying it by 2irr Jo(k±r) exp(iwi) and integrating over r and t, 
one founds 



where Eqs. apply, and 



~ ^j,C = X,[P\{a,P)S^"^-^', (39) 



X/[P](a,p) = kaB[P](a+^,Vkk"a). (40) 



k 

Eq. H39fl can be readily integrated with the boundary condition C = at z 



with 



C(a,[3,z) = ^XAP^-^le^^-^-^, (41) 
9 



^ a + ^-^ + W ] - (42) 



FIG. 2: Sketch of the generated spatio-temporal spectrum (symmetrized for k± < 0) when (3 > 0. The straight lines are slash 
and backslash spectra. 



Eq. (|41|l can be interpreted in the (u>, k±), or (a, f3), planes and shows that, for large propagation distances, C tends 
to a Dirac 6 centered at g = 0: the propagation acts as a spatio-temporal filter, selecting specific combinations of 
frequencies and wavevectors. 

The condition g — 0, in the {uj, kj_) plane, is 

Y {uj -^ ] -2k=W (43) 

which gives a hyperbola, as sketched in figure El It is apparent that Eq. (|43J) is the dispersion relation corresponding 
to Eq. ®, therefore the resulting pulsed beam is a progressive undistorted wave traveling at inverse differential 
velocity [3. In reported experimental results there is a clear evidence of this spectral hyperbola, see e.g. as well 
as of a structured spectrum. |10ll47| as that corresponding to the splitting described in the following. 

The asymptotes over which energy is concentrated correspond to slash and backslash X-waves. Thus, in any 
nonlinear process which can be reduced to equation l|35|) . X-wave packets are spontaneously generated. I like to stress 
that one of the arms always starts at u> = 0, while in the case f3 — Eq. (|43|l yields the exact X-shaped spectrum. 
The specific features will in general depend on the source spectrum , as elucidated below with an example. 

Note also that the same analysis holds true, with simple changes, when \Pnl(z, t, r) — V(t — /3z, r) exp{ik z z) 1 with 
k z depending on the specific nonlinear process (e.g. k z is the wavevector mismatch for optical harmonics generation), 
detailed cases will be reported elsewhere. 



Energy of the generated X-wave 



The energy distribution function is, after (|22|L 

f0 ° 4TT 2 \ C(a,P,z)\ 2 r° 4^ 2 |A/[P]| 2 (a,/3) sinjg z^ 2 

ka J Q ka g 



W)= l - 1 t„ da = I T da - ( 44 ) 



As z — > oo 

sin(gz) 2 2ir z 



9 2 
with 

__ /?(/?- 2£) 



5(a-a), (45) 



2k" 09 - (3) • 

Hence for large z 



(46) 



Rw^kn 7 ft 

^ ^p\W P ^+ | 7 '^ 77 «)| 20 o(«), (47) 
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with 6*o the unit step function. 

The generated X-wave packet grows with a linear efficiency with respect to the propagation distance. In the case 
(3 7^ 0, from the condition a > in (|4*7I) it appears that energy is distributed in the intervals (3 < and (3 < (3 < 2/3 
when /3 > 0, and in (3 < (3 < and (3 < 2(3 when (3 < 0. This implies that £p does not have a continuous support (see 
the example in figure|3| discussed below) For large |/3|, with sufficiently separated domains along (3, this is expected to 
provide pulse-splitting. In fact, some components of the generated X-wave will travel at a velocity sensibly different 
from the pump wave, and after some propagation length satellite packets may appear (such as, e.g., those discussed 
in jH). 

There are two relevant limits for Eq. I|47|) . 

As (3 — > 0, i.e. the pump travels at the linear group velocity of the medium, and a — > —f3/2k", it is 

£^^\B[P](^,-Vw^)\ 2 e (-(3). (48) 

Eq. I)48|) shows that only the portion of the pump spectrum laying on the X-spectrum, i.e. on the line k± = —\/kk"uj, 
contributes to the energy of the generated beam. The propagation filters out all of the spectral components that do 
not belong to an X-wave. 

The other relevant limit corresponds to a pump velocity largely different from the linear group velocity of the 
medium: — * oo. From a physical point of view, the nonlinear interaction is strongly hampered by the generated 
beam which rapidly gets separated from the pump, producing a typical "Cerenkov emission" halo. From 147|) it is 

(a - -f3/k") 

* 8 -^W lB[m -^F)| 2 *o(-/?)- (49) 

For large \/3\ the X-wave generation is strongly inhibited, and exclusively the on-axis spatio-temporal spectrum plays 
a role. Only the branch (3 < is relevant in this limit. Furthermore, because of the factor |/3|, the energy distribution 
is expected to be peaked at some value (3 < 0, as shown below. 



Example: Gaussian pump beam 

It is instructive to examine an example, by assuming a Gaussian lineshape for the source: 

B[P}(k U w)=B$exp(-^--^), (50) 

such that Bq measures the pump beam fluence, and Awp and Akp provide the on axis temporal and spatial spectra, 
respectively. 

In the case [3 = (source traveling at the linear group velocity at luq) the energy distribution function, after Eq. 
63, is 

4ir 3 B?:kz , (3 2 , „. 
Sfi - fc ° exp(--^)g (-/3) (51) 

showing that the generated spectrum encompasses slash X- waves faster than the source {(3 < 0) , and spread in velocity 
with the characteristic value 

A ^ p = [ 4(P)^| + 4F^4 ]_1 ' (52) 

such that the more the source is spectrally (in space or time) narrow the more "ideal" (i.e. un-spreading, A(3p small) 
is the generated X-wave-packet. 

In the case (3 ^ the situation is more complicated. From (|47|l and l|50l) . it is £p = when a((3) < 0; while for 
other values of (3 



A^Blkz(3 2 -2(3j3 (3 2 (3 2 - g P (f3P - ^ 2 ) 1 

E < 3 = w> w^w eM ~w^w — m — ] (53) 
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FIG. 3: Normalized energy distribution W Vs inverse differential velocity (3 for two different values of pump velocity determined 
by p. (a) branch /3 < for p = 0.01 A/3 P ; (b) branch p < p < 2Ap P for p = 0.01A/3 P ; (c) as in (a) with p = 0.5A/3 P ; (d) as 
in (b) with P = 0.5A/3p. All quantities are dimension-less. 



with q P = kA{3p/k"Ak P . 

In figure 01 1 show an example of the application of Eq. I|58[l. by reporting the normalized energy distribution 
W — £pk" /A-K^B^kz. The pump beam, with uiq such that the wavelength is Ao = 1/im, is chosen with Awp = 2tt/ At, 
with At = 500/s, and Akp = Adko with AO — O.lrad. For the medium parameters no = 2, k' = c/no, k" = 
360 x 10~ 28 s 2 /m. As a result, p = 10~ 12 s/ra, with a velocity bandwidth = A/3 P = 10 4 m/s, and g p = 10~ 4 . 

In figureOl ( a ) an d (b), the two branches are displayed for /3 = 0.0lA/3p, corresponding to a pump beam velocity 
of Vp = 0.499999c. In this case, the energy distribution tends to that at = 0, and no pulse splitting is expected 
(the two branches are almost contiguous). 

In figure |3 (c) and (d), the two branches are shown for (3 = 0.5A/3p, corresponding to a pump beam velocity of 
Vp = 0.499966c. Two distinct lobes are present, with the most energetic peaked around = —2(3. In this case a 
satellite X-wave packet is expected to separate from the pump beam in propagation. 

The X-wave will in general split for sufficiently long propagation, following the peaks of £p. At higher orders 
the resulting packets will become new sources for the perturbation equation, and additional splitting will appear, 
consistently with the numerical results (see e.g. |35H45) l. 

This is expected in second- harmonic generation, initially considered in |l5||. where /3 is the temporal walk-off between 
the two harmonics (details will be given elsewhere); as well as in dispe"ve propagation, as described in 0, and 
investigated experimentally in |4f| . The spectral widening can indeed provide sources traveling at a group velocity 
different from that at wo 

In figure |21 note the difference in the vertical scale between panels (a), (b) and (c), (d). In all cases the pump 
velocity is about 1/k', showing that the mismatch between the pump and the linear group velocities at ujq drastically 
affects the generated X-wave spectrum. 

With regards to the efficiency of X-wave generation, in the case = 0, by integrating Ij48(l it is found 

£ = 27TV [ B } kZ Ap P = 4^Bg k z , ^ P . (54) 
k" I . y k Aul 
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The energy grows with A/3p, i.e. with an extension of the spatio-temporal spectrum of the pump. After Ea. ll54() 
one founds that, for a wide spatial spectrum (i.e. large Akp, corresponding to a tightly focused pump beam), the 
efficiency is mainly determined by the on- axis temporal spectrum. Conversely, for large on-axis spectrum Atop, the 
energy of the generated X-wave grows with Akp. 
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The other relevant limit is (3 » Af3p, such that the resulting energy distribution is solely determined by the 
branch (3 < 0, and it is 

£= 2^ _ z Ak P 
P 

Thus the efficiency is mainly affected by the spatial profile of the pump beam and it goes like 1/|/3|, consistently with 
Eq. Pty. 



HIGHLY NONLINEAR REGIME: AN EFFECTIVE ID NLS 



Even when using some working hypotheses, it is possible to deepen the previous study and consider the highly 
nonlinear regime. I will consider the multidimensional Schroedinger equation: 

1 k" kn^ 

id z A + ik'd T A + —V 2 A - —d TT A + —^\A\ 2 A = 0. (56) 
Ik " I hq 

This ma y p rovide the correct trend observed in the experiments, including quadratic nonlinearity in certain regimes, 
(see e.g. j53) as long as envelope shocks (see |45jV or higher order phenomena such as plasma formation play a 
negligible role. 

Writing A as a X-wave expansion (0 < a < 00 and —00 < [3 < 00): 

A = J X(a,f3,z)J (Vkk T ar)e^ a+ ^ k " )(t - l3z) dadp (57) 

and taking the Fourier Bessel transform of Eq. I|56|) . evaluated at k± = \J kk"a, and u> = a + /3/k", the coupled 
X-wave equations are derived: 

id z X(a, (3, z) - £rX(a, (3, z) + ^Q{a, (3, z) = 0. (58) 

The nonlinear polarization is 



Q = X/[\A\ 2 A*}e- l{a+ ^ )l3z = 

kaj K.(a,d)X(ai,(3i,z)X(a2,(32,z)X(a 3 ,P3,z)* (59) 
5[a + a 3 - ai - a 2 + {fi + 03 - Pi - fo)/k"]B(a, (3, a, (3, z)dad/3, 
with a = (a 1 ,a 2 ,a 3 ), (3 = {(3i, P2, fts),]^ 

/•oo 

/C(a, ai, a2, a 3 ) — / Jo(V kk"ar) Jo(V kk"a\r) Jo(V kk" a 2 r) Jo(V kk" a 3 r)r dr (60) 
Jo 

and 



9 = exp{i[-(a + f3/k")(3 + + (3 x jk")(3 x + (a 2 + (3 2 /k")f3 2 - {a 3 + p 3 /k")(3 3 ]z} . (61) 

A remarkable result is obtained if the solution of l|58|l is approximated by an X-wave-packet which, in order to 
avoid unnecessary complexities, is centered around the medium group velocity, i.e. (3 = 0. This is standard in usual 
coupled mode theory; [59j in the present case a "mode" is an X-shaped 3D wavepackct. 

As discussed below, if z is sufficiently small, it is possible to take = 1 in l|59() : in this regime I write X(a,(3) = 
f(a)C((3) in order to represent the nonlinear modulation (with envelope C) of an X-wave. For definiteness, I consider 
/ = f p , where f P (a) is the spectrum of a basis X-wave. If C is peaked around (3 = 0, all components travel nearly at 
the same velocity and = 1. Multiplying by At: 2 f p (a)/ka and integrating over a, from Eq. (|58|l one obtains: 



td z C(f3, z) - £rC{J3, z) + ^j x ((3 + p 1 -(3 2 - /3 3 )C(A) C{(3 2 ) C{f3 3 )*dp = 0. (62) 
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with the interaction kernel: 

X(j) = An 2 J fC(a, ai, a 2 , a 3 )f(a)f(ai) f(a 2 )* f{a 3 )*S[a + a x - a 2 - a 3 + jj;]da. (63) 

Eq. H62fl is a Zakharov equation, and taking the Fourier transform of C (see (G2J)) one obtains the ID nonlinear 
Schroedinger equation (NLS): 

d c k" d 2 c kn 2 , . , l9 ,„ 
oz 2 o s z n 

Hence the evolution of an X-wave packet in a nonlinear Kerr medium can be approximated by an effective 1+1D 
nonlinear Schroedinger equation with a non-homogeneous nonlinearity. The latter, given by a(s) with dimensions of 
an inverse area, is expressed by the Fourier transform of the interaction kernel \ : 

/oo 
Xh)e ns/k " 'drf. (65) 
-oo 

After some algebra it can be written as 

p CO ^ 

o-(s)= |27iVfc 7 V / (p) | 4 rdr, (66) 
Jo 

which is the spatial overlap of the component X-wave profile at — 0, with 



V/ (p) =(p/ (p)(r ) g)= f p (a)J (Vkk"ar)e~ las da. (67) 
Jo 

Taking for c a solution of Ipjlfl. and C its Fourier transform according to l(3*2l . A reads 

A= [ C{0,z)ip, ( - p \r,t- 0z,0)d0 = [ c(s, z)£/ (p) (r,i, s, z)ds, (68) 



£/ p) (r,<,s,z) = -L= f e - l[ - s+t) ^ +l & z V/ ^\r,t- 0z)d0. (69) 



with 

£,V>(r.t.a.z) = - 

Ik" 

can be considered the basis for the wavelet transform of A with respect to the on-axis temporal variable. 60| 
If 8 = 1 the spreading due to the different velocities of the component X- waves (given by the quadratic terms in in 
the exponential in l|69l) ) is negligible, and £ can be approximated by its value at z = 0: 

^(r, t, s, z) S t, s, 0) = 2irVk*5(t + s)^(r, t) (70) 

and l (58j l becomes 

A = 2irv / ~k"c(-t,z)ip / ( p) (r,t). (71) 

The 3D+1 nonlinear evolution problem H56J) reduces, under suitable approximations, to a 1+1D model l|M[l . In the 
following I will address the hypotheses leading to this result and to its straightforward consequences. 

CONDITIONS FOR OBSERVING SOLITONS 

Eq. is valid as far as 6 = 1 in JEjJ) . Writing 

Q = exp(i6z), (72) 

with 

= ( ai + 1 /k")0 1 + (a 2 + [3 2 /k")f3 2 ~(a + 0/k")0 - (a 3 + 3 /k")0 3 , (73) 
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it must be \9\z « 1. 

9 is a quadratic function of /3s; if 5{3 > is the velocity bandwidth ((3 < \5j3\), and observing that the spectrum of 
the fundamental X-waves decays as exp(-aA) (thus, roughly a < 1/A), the validity of this approximation leads to 
a constrained maximization of 9. After a simple analysis, it must be 

z « TifliT = z x (74) 

max(|t/|) 

while, if 1/A < 25/3/k", the soliton regime: 



2A- 



and, if 1/A > 25/3/k", the chirp regime: 



45/3 
~A~' 



(76) 



zx 

As expected, the smaller 5(3, the larger zx- 

The above can be recast in a much simpler and more insightful formulation by introducing the dispersion (Ldisp) 
and the diffraction (Ldiff) lengths. 

If 7o = k" /25(3 is the s-width of c, which in general does not correspond to the on-axis temporal duration of the 
pulsed beam (this holds approximately only in the soliton regime, see dJ), it is Ldisp — /k" . 

If Wo = 1/Afcj_ is the beam waist at the center of the pulse, with Ak± the spatial spectrum, by the general 
properties of X-waves it is (a < 1/A) Wo = A/Vkk" , and Ldiff — fc W 2 - 

When A is much smaller than Tq, the X-wave dominates the whole spatio-temporal beam profile, this is the chirp 
regime. Conversely, when Tq < A, the soliton regime is attained. In this case, the on-axis temporal duration is also 
affected by the envelope c, while the X-wave mainly determines the spatial profile (see l|71|l). 

For instance, assuming to be in the soliton regime (as witnessed by non-trivial nonlinear dynamics), taking a 
beam waist of Wo = 50/i m and an on-axis temporal duration Tq = 200/s it is zx — 8cm, taking uq = 2, k" = 
360 x 10 _28 s 2 /m and Ao = 500nm, a value far beyond the propagation distances reported in experiments. [7|,|46j zx 
is typically of the order of the smallest between Ldisp and Ldiff, because it measures the distance after which the 
X-wave packet starts to spread due to its finite energy. Therefore one expects the approximations leading to the ID 
NLS model Ea. (|64|l to remain valid at least in the early stages of the dynamics, a few centimeters in experiments. 
When z increases, the nonlinear response gets reduced because of the sliding between the component X-waves, and 
the propagation becomes in essence the linear evolution of the generated X-wave pattern. 

Once the validity of Q64JI is ascertained, the corresponding nonlinear dynamics is determined by considering the 
the effective nonlinear length Lxl — n-o/(kn2 c(0) Cq), with Co = |c(s — 0, z — 0)| 2 , which relates to the input peak 
intensity, omitting inessential numerical factors, by c\ = Tq W 2 (i.e. c\ is approximately given by the peak power, see 
below). In the definition of Lnl, the value of cr(s) is considered at s = 0, because prominent nonlinear effects are 
expected near the peak of the pulsed beam. 

The depth of focus of the progressive undistorted wave is such that the on-axis intensity keeps constant in propa- 
gation, thus enhancing nonlinear effects. This is quite similar to the nonlinear dynamics of guided modes, where the 
tightly confined light may favor nonlinear phenomena. 

When the temporal duration Tq of the envelope c is smaller than the duration of the modulated X-wave, approxi- 
mately A, as far as Ldisp >> Lxl the nonlinearity will mainly play as a self phase modulation. When Ldisp — Lnl 
the so called N = 1 fundamental soliton solution (with a sech profile for c, and N the number of eigenvalues in the 
inverse scattering problem of H64J) . with <r(s) — > cr(0)) is attained. However, it can be readily seen that Ldisp > zx 
(in the previous example Ldisp — lwz), hence the X-wave dynamics dominates the dispersion. Even if a N = 1 soliton 
emerges, it essentially behaves as a non-dispersing non-diffracting X-wave. 

The situation can be completely different when the peak power increases and Lxl reduces. In this case the 
parameter N, given by y/ Ldisp I Lnl, becomes greater than unity, and hence higher-order solitons and breathers are 
foreseen. They will be discussed below. 



FUNDAMENTAL X- WAVE/FUNDAMENTAL SOLITON 

Under suitable conditions, satisfied in the early stages of propagation, the nonlinear dynamics of an X-wave in a 
Kerr medium can be described by the ID nonlinear Schroedinger equation with a non-uniform nonlinearity profile 
<r(s), determined by the spatial (self-)overlap of the progressive undistorted wave. 
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FIG. 4: Nonlinearity profile induced by the spatial overlap of the fundamental X-wave. 



As I have shown earlier, the whole X-wave space is spanned by the fundamental X-waves. The fundamental 
X-shaped profile, integrable with respect to r and typically observed in simulations and experiments, is given by 
fo(a) — (Vk/n) Aa exp(-Aa) (see l |27t ). I will consider the effective ID NLS with such spectrum. It can be found 
that 



\<P/ 



(0)|4 _ 



fc 2 A 4 



(s 2 + A 2 



[(s 2 + A 2 ) 2 + (jfe k" r 2 ) 2 - 2(s 2 - A 2 )k k" r 2 ] 3 ' 



which gives 



a(s) 



8k" k s_ 

,2 a >■ 



(77) 



(78) 



5A 2 Juy A- 

f a , plotted in figure 01 and such that / CT (0) = 1, can be evaluated either numerically or analytically. In the latter 
case it can be expressed as [61| 



fM = + 3i(l + M 2 )>g(~ £) - log(&) - log(^l!) + log(-i^)]} 



(i+^ 2 ) 



2/i ' 



-2/i' 



2/i 



2/i 



(79) 



Therefore, a(s) is a bell-shaped function. Observing that its spatial extension is of the order of A, when the 
on-axis temporal dynamics is dominated by the c envelope (i.e. the temporal width of c is smaller that A), it can be 
approximated by u(0): 



a(0) 



8 k k" 
5A 2 



1 



The last relation stems from the fact that the spatial spectrum is Ak^_ 
thus be approximately rewritten as 



(80) 

1/Wo = \/fcfc 77 /A. The effective NLS can 



.dc k" d 2 c kri2 
''d~z + + 7i W 2 



! c = 0, (81) 



so that the effective nonlinear coefficient has an additional factor, roughly given by 1/Wq (or equivalently Ak 2 ), if 
compared to the plane wave value. This shows the compensation of diffraction owing to the X-waves, which behave 
as "modes of free space". c(s) resembles the mode-overlap in a waveguide and the corresponding support to the 
nonlinear response along propagation. This holds valid as long as z << zx, i-e. as long as the nonlinearity is not 
averaged out by the "sliding" of the X-waves. 

The fundamental soliton (when <r(s) = er(0)) can be expressed in terms of the peak Cq of the energy distribution 
function; it reads 



.4 



l c 2 kn 2 a(0) 



.Cq fcn 2 cr(0) 



c sech(^/ ~" '\'J t Z s ) exp[ C" " K "' z]C/ (0) (r, t, z, s)ds, 



k" no 



2n 



(82) 
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which, using (|71|l and i|80|) . can be approximated by 

A = 2 7 rVF Co sech(,/3^i) exp(i^^z)^ (r, t). (83) 
y Wo ft" n 2noWo 

As anticipated, from the expression of if/^\ the peak intensity 1$ relates to c$ by Cq = 2oWq (omitting inessential 
factors). 

Eqs. H82fl and l|83[l show the "dressing" mechanism associated with the nonlinearity: the latter acts only on the shape 
of the envelope which, even in the linear limit, would travel almost undistorted (as far as z « zx). Since this analysis 
can be repeated for any X-wave, nonlinear X-waves are not numerable families of solutions (like multidimensional 
solitary waves), but have the same "cardinality" of the "space" of linear X-waves. 

Conversely, in spite of this remarkable difference with solitary waves, the (approximate) validity of an integrable 
model such as l|81fl seems to establish a strong link with solitons. Clearly this analysis is far from being a rigorous 
settlement of this conclusion, but could stimulate further research. In the following section I will discuss some 
additional consequences of (|64|l and l|81|) . pointing out non-trivial nonlinear dynamics of X-waves. 



NONLINEARLY CHIRPED X-WAVES AND COMPRESSION 



The fundamental X- wave/fundamental soliton discussed above is of limited physical interest, because its invariant 
propagation is already "embedded" in its linear counterpart. If C{f3) is a very narrow function, or equivalently the 
s— width 7o of c(s, z) is much wider than A, the on-axis profile of the 3D-beam is mainly determined by the X-wave. 
This has been indicates as the chirp regime. 

In Eq. lt)4|) the second derivative can be neglected, and the solution is essentially self-phase modulation: 

c(s, z) = A eiq)[i—?-cr(s)Alz}. (84) 
Proceeding as before l|84|l yields the approximate nonlinearly chirped basis X-wave 

A - 2 7 rVfc A /|^ (p) M) ex P [z^ f a (t)). (85) 
V Wo n l 

There is a fundamental difference between this type of chirp and what is typically considered in fiber propagation, 
|62l | determined by the temporal power profile |c| 2 . Even an idealized X-wave (corresponding to C((3) proportional to 
a Dirac delta with respect to (3, and hence to a constant c) is chirped in a nonlinear medium because of its spatial 
profile ( reflected in f a (s)). 

Since the earliest work on pulse splitting in normally dispersive Kerr media, it has been known that, at the initial 
stages, the on-axis temporal spectrum exhibits a double peak. |37j Figure shows the instantaneous frequency 

kn 2 z df a 

5w = ^T-TT- (86) 

uqIq at 

For higher order basis X-waves, more complicated spectral modulations can be expected. It is natural to identify 
the origin of this spectral splitting as a consequence of the self-phase modulation considered here. It is also well 
stated that a chirped pulse may compress when propagating under anomalous dispersion. With reference to the two 
frequencies at the peaks in figure they propagate according to Eq. (|64[) with opp osite velocities, such that the 
pulse gets compressed in the temporal domain. In the pulse splitting considered in |37| . on-axis temporal compression 
immediately follows the spectral splitting. 

Thus, taking the whole process of pulse-splitting and admitting that an X-wave is formed in the initial stage (as in 
[13) i the origin of the spectral splitting can be attributed to self-phase modulation of an X-wave. Nonlinear processes 
are therefore well suited for generating chirped X-waves, such as those recently addressed in 



SPLITTING AND REPLENISHMENT IN KERR MEDIA AS A HIGHER ORDER SOLITON 

While the fundamental soliton and the chirp provide the simplest corollaries of Eq. I|64|) , the self-trapped behavior 
of the 3D beam is fundamentally due to the X-shape. More interesting dynamics can be described referring to multi- 
soliton solutions (obviously in the soliton regime). The N > 1 solitons [62| are natural concepts in explaining splitting 
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FIG. 5: Instantaneous frequency corresponding to X-wave self-overlap in Fig. 0] 

and replenishment, investigated numerically in |35| and experimentally in [jij ]. Noteworthy, breathing linear X-waves 
have been reported. |6J, l65| 

Even if the replenishment dynamics has been examined with reference to water while including additional terms 
to Eq. H56JI . the most relevant features are somehow taken into account by this "simple" model, with plasma forma- 
tion and higher order dispersion playing a perturbative role. 35] Highly nonlinear phenomena, such as the shocks 
considered in |45|. are expected at very high fluences and will be neglected hereby. 

For example, the N=2 soliton j^] provides the approximate breathing nonlinear X-wave 



k" 

A = 2n= 



no 



t \ < ,^-r-'^/ (0) ^ 

T V kn 2 a(0) L D T 



2fc"W fm z t A(t + iA) 

^Tv~ V ~2 U[ L^' V [(t + iAf - k k" r 2 ]3/ 2 ' 



(87) 



with 



^ cosh(3r) + 3 exp(4i£) cosh(r) 
cosh(4r) + 4 cosh(2r) + 3 cos(4£) ' 



it/2 



Figure |B| shows a typical spatio-temporal profile obtained after Eq. 1|87|) . The periodic depletion and replenishment 
of the X-shaped distribution is apparent. In propagation, the two-solitons, or breather-, solution pulsates and the 
beam evolves retaining most of its energy localized, but exhibiting the non-trivial nonlinear dynamics of the X-wave. 
The generation of a breather can also be thought at the origin of splitting in the first stages of propagation. It is also 
noticeable that higher order solitons exhibit the spectral splitting typically described in numerical simulations. 

Before concluding, it is fruitful to summarize the picture of the splitting/replenishment phenomenon. At the 
beginning a wide bell-shaped pulsed beam evolves into an X-wave, owing to the spatio-temporal pattern formation of 
X-wave instability 01 . Then, an effective anomalous dispersion is experienced by the envelope of the finite-energy 
X-wave packet, and the chirp regime is entered with spectral splitting and on-axis compression. Once the envelope 
width is sufficiently reduced, the soliton regime starts, and the increased intensity through compression feeds the 
generation of a higher order soliton, or breather. After some spatio-temporal oscillations, several mechanisms may 
intervene to stop the periodic behavior, e.g. losses (eventually of nonlinear origin, such as two-photon absorption) or 
simply that, for large propagation distances, the nonlinear response average out due to the sliding between components 
of the finite energy X-wave packet. 



CONCLUSIONS 



Progressive undistorted waves, a generalization of the non-monochromatic realm of self-invariant beams, seem to 
emerge as a valuable tool in numerous fields of applied research, from telecommunications to biophysics (supported 
by the recent experiments in water). Their natural appearance in nonlinear processes, analyzed hereby in the optical 
domain but also available in acoustics as well as in Bose Einstein condensation, can be considered a fundamental 
result. The formation of X-waves during frequency generation appears as natural as the use of their paradigm in 
interpreting 3D+1 nonlinear dynamics. 
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transversal coordinate x 
(arb. units) 



FIG. 6: (Color online) Typical spatio-temporal profile (at y = 0) of a breathing X-wave. Two isosurfaces are displayed: the 
darkest (red) corresponds to higher intensity. 

While X-waves are known to carry infinite energy, and their superposition has been previously employed to build 
finite energy solutions, in this paper for the first time I have attempted to establish a general -albeit preliminary- 
theoretical framework for the extension of this approach to the nonlinear case. 

The various steps of a basic nonlinear process, such as pulse-splitting in normally dispersive media, can be revisited 
in terms of X-wave nonlinear dynamics: from the initial pulse compression to splitting/replenishment. The results 
open the way to the investigation of elastic collisions between nonlinear X-waves, as well as to establishing a link with 
parametric solitary waves in quadratic media. 

According to standard textbooks, a soliton (or a bullet) is a non-perturbative solution of a nonlinear wave equa- 
tion. In this sense, the non dispersive and non diffracting wave-packets analyzed here cannot be considered solitons, 
because they exists even when a nonlinear susceptibility lacks. Nevertheless, an intriguing connection between pro- 
gressive undistorted waves and solitons seems to be at the inner stem of numerically and experimentally investigated 
phenomena. In this respect nonlinear X-waves can be considered a sort of chimera of modern nonlinear physics. 
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